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Fractal surfaces (’patchwork quilts’) are shown to arise under most general circumstances involving 
simple bitwise operations between real numbers. A theory is presented for all deterministic bitwise 
operations on a finite alphabet. It is shown that these models give rise to a roughness exponent 
H that shapes the resulting spatial patterns, larger values of the exponent leading to coarser surfaces. 
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I. INTRODUCTION 


Fractal surfaces mi are ubiquitously found in biological and physical systems at all scales, ranging 
from atoms to galaxies Mathematically, fractals arise from iterated function systems mi, 

strange attractors [1 , critical phenomena uni, cellular automata [nUU, substitution systems mm 
and any context where some hierarchical structure is present mm- Because of these connections, 
fractals are also important in some recent approaches to nonequilibrium statistical mechanics of steady 
states (T6HT9] . 

Fractality is intimately connected to power laws, real-valued dimensions and exponential growth of 
details as resulting from an increase in the resolution of a geometric object [20]. Although the rigorous 
definition of a fractal requires that the Hausdorff-Besicovitch dimension Dp is strictly greater than the 
Euclidean dimension D of the space in which the fractal object is embedded, there are cases in which 
surfaces are sufficiently broken at all length scales so as to deserve being named ‘fractals ’ mcaaEu. 
Most of these surfaces are random and governed by growth, adsorption and deposition processes on 
interfaces [20] 1221428] . 

The multidisciplinary field of disorderly surface growth has experienced a rapid development [23]. 
However, fractal surfaces formed by deterministic processes and hierarchical rules are also interesting 
and can be useful, e.g., as models for cityscapes [20]. In this article, by regarding ‘fractals’ in the 
broad sense mentioned in the previous paragraph, we construct a wide variety of such surfaces and 
show that they possess fractal self-affine properties. If a surface F(x , y) is self-affine it satisfies |22j [23] 

F{x, y) ~ b~ H F(bx, by) (1) 

where H G M is the roughness exponent characterizing the self-affine scaling [TJ [29] . We prove that the 
surfaces here constructed obey Eq. 0. Although our construction proceeds abstractly, we illustrate it 
with specific numerical examples, and we believe that the generality of the method is such that it may 
find many applications in the modeling of complex physical systems. Although the notation may seem 
unfamiliar, the mathematics behind is elementary, and is based on generalized bitwise arithmetic on 
real numbers. To introduce this idea, let us consider two real numbers a and b (that we may truncate 
to a finite number of digits after the decimal point). If we expand a and b in base 10, we see that the 
ordinary sum of these numbers splits into the ordinary sum of two different parts 

cl + b = (a +io b) + (a+io&) (2) 

here +io denotes addition modulo 10 and +io denotes the contribution of the carries to the sum a + b. 
If 0 < a < 10 and 0 < b < 10 then, it is clear that a + b = a +io b if a + b < 10. If a + b > 10 then 
cl +io b = a + b — 10. If a and/or b are larger than 10, the sum a + b is clearly reduced to separately 
considering the digits of a and b and adding them, taking care of the carries. Addition modulo 10 means 
that the carries are neglected. For example, if a = 5.6782 and b = 3.6754 we have that a + b = 9.3536. 
This sum can be seen as the ordinary addition of two different contributions a +io b = 8.2436 and 
a+iofr = 1-1100. The structure of the operators +io and +io is very interesting. All positions within 
the numbers are independent of each other as regards the bitwise action of +io: Each two digits 
corresponding to the same power of ten are added modulo ten and no carry is transferred from one 
position to another. Therefore, under the action of such operator, the positions within the number are 
‘uncoupled’ and independent. Since each position within a number in a standard positional number 
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FIG. 1: The functions f(x, y) = x + y (left), g(x, y) = x +2 y (center) and h(x , y) = x+ 2 V (right) for x £ [0,1] and y £ [0,1]. We 
have f(x , y) = g(x, y) + h(x , y) everywhere. Dark blue means a zero value, green a value of 1 and dark red a value of 2. 


system corresponds to a different power of the base, adding modulo 10 two real numbers means 
performing the same operation at all scales, the latter being independent of each other. We claim 
that, as a result of this, the function g(x, y) = x +10 y (where x and y are real numbers) has fractal 
features. It is a discontinuous and non-differentiable function that exhibits the same details at all 
scales, displaying self-similarity. In fact this is generally the case for the function g(x, y) = x + p y, 
where p 6 N, p > 2 is any base, with + p denoting addition modulo p MED). We shall prove 
below that such a function g(x, y) obeys Eq. 0' Furthermore, since f(x,y) = x + y is not a fractal, 
the function h(x, y) = x+ p y = x + y — (x + p y) must be a fractal as well so as to compensate the 
discontinuities of g(x,y). In Fig. 0 the functions f(x, y) = x + y (left), g(x, y) = x +2 y (right) and 
h(x, y) = x+ 2 y are plotted for x € [0,1] and y £ [0,1]. While f(x, y) = x + y is just an inclined plane, 
the functions g(x, y) = x +2 y and h(x, y) = x+ 2 y exhibit a more complex behavior displaying ‘squares 
within squares’ at all scales. The behavior of these functions is such that under ordinary addition, all 
discontinuous jumps disappear at all scales, yielding the smooth function f(x, y) (which is everywhere 
equal to g(x, y) + h{x , y)). 

The purpose of this paper is to construct all possible bitwise operators acting on a finite alphabet of 
p symbols (to which + p and + p do belong, and which we explicitly construct as well, as an example) 
and to elucidate their associated fractal features. We thus find a new kind of surfaces with self-affine 
properties, that we term ‘patchwork quilts’. We constructively prove the general fact that bitwise 
operators locally acting on each point (x, y) in the plane M 2 leads to these surfaces, providing explicit 
mathematical expressions for them. The results can be easily extended to hypersurfaces in M n . We 
explicitly construct the generalized bitwise operators involved in this process by making use of the 
framework of digital calculus, that we have very recently introduced H2IE2M. We show that these 
abstract surfaces give rise to a roughness exponent H that shapes the resulting spatial patterns, larger 
values of the exponent leading to coarser surfaces. 

Recently, we have found a new construction for fractals based on a fractal decomposition of a given 
function [33] linking the resulting fractal objects to finite group theory [33] . The construction of fractal 
surfaces given here is new and different to our previous one in [33]. Although operators acting on two 
variables are again considered, these operators are not restricted to have the latin square property as 
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and the construction is different even when the concepts that are used in it are closely related. 

In Section [TT] we introduce the main concept of digital calculus, the digit function [32, [33j and we 
construct all discrete operators acting on two variables whose values are restricted to be any of the p 
integers in the interval [0,p— 1]. With the aid of these tools, we then construct local bitwise operators 
acting on points in and we rigorously prove that all give rise to (hyper)surfaces satisfying the 
self-affine property meant by Eq. 0- We then also prove that the roughness exponent behaves as in 
experimental physical systems, larger values of it leading to coarser surfaces [29, |3U GEJ. hi Section 
TIT we illustrate these results with specific numerical examples. 


II. DIGITAL CALCULUS: THE DIGIT FUNCTION, DISCRETE OPERATORS AND 
GENERALIZED BITWISE ARITHMETIC 


In this article S denotes the set of p integers in the interval [0,p — 1]. The digit function d p (k,x) 
for p 6 N, k e Z and x G M, is the surjective mapping M -» S defined as na Eli E3] 


d p (k,x) = 


.P 


p 


p 


x 

k+1 


(3) 


and gives the k-th digit of the real number x (when it is non-negative) in a positional numeral system 
in radix p > 1. If p = 1 the digit function satisfies di(fc,x) = di(0, x) = 0 and it does not relate to 
a positional numeral system. In Eq. ([3]) |_.. .J denotes the floor function (lower closest integer) of the 
quantity between the brackets. 

With the digit function, we can express any real number x as P31E21E3] 

U°gp Nil 

x = sign(x) E p k d p (k, |x|) 

k=— oo 

The integer part of x is given by 

U°gp MJ 

sign(x)[|x[J = sign(x) ^ p k d p (k,\x\) 

k =o 


(4) 


(5) 


Example: In the decimal radix p = 10, the number it = 3.1415 ... has digits dio(0,7r) = 3, dio(—1, n) = 
1, di 0 (—2,7 t) = 4, di 0 (—3,7r) = 1, di 0 (-4,7r) = 5, etc. 

In general, a truncation to D digits after the radix point is given by 


sign(x)p = sign(x) 


U°gp MJ 

E J 

k=—D 


p k d p (k,\x\) 


( 6 ) 


We shall call p~ D \p D .. .J, with D being an integer number, the coarse graining operator. 

Most properties of the digit function are easily understood from its plot, as shown in Fig. [2] for the 
case p = 3 and k = —2, —3, —4. The digit function dp(k, x) is in all cases a periodic staircase that is 
set to zero each time that x is an integer multiple of p k+l . The digit function possesses an important 
scaling property [33] that we shall use below. We have 


d p (k,p m x) = 


p x 

pk 


p 


p x 

qr\k- 1-1 


X 


,k—m 


P 


P 


X 


k — ITL+l 


P 


= d p (k — m, x) 


(7) 
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FIG. 2: The digit function d 3 (k,x) plotted for k = —2, —3, —4 (from top to bottom) in the interval x E [0,1]. m 


Let n be a natural number. Then, the following relationship also holds 


d np(k,x) = dp (fc,T) +pd n 


= d « ( ^k) + nd P (*> 


( 8 ) 


This can be easily seen by using the definition, since we have 


n" 


d p[ k ’-r) + P d n k, 


( /c+l 


P 


p k n k 


■P 


pk+l n k 


+ P 


X 


pk+l n k 


np 


pk-\- \ j^k-\-\ 


( np) k ^ ( np) k+l 


= d np (k,x) 


(9) 
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which proves Eq. Q. 

Let us now construct the operators + p and + p mentioned in the introduction. We consider a and b 
nonnegative real numbers, for simplicity. From the digit function, Eq. Q, we have that, on one hand 

U°gp(“+&)J 

a + b= 22 p k d p (k,a + b) ( 10 ) 

k=— oo 

and, on the other hand, 

Liogp a| L lo s P b \ max ( L lo g P a \ , L'ogp b \} 

a + b= 22 P k d p (k,a)+ 22 p k d p (k,b) = 22 p k {d p (k, a) + d p (k,b)) (11) 

k=— oo k=— oo k=— oo 

Since p 2 > d p (k,a) + d p (k,b), we have d p 2 (0,d p (fc, a) + d p (k,b)) = d p (k,a) + d p (k,b) and we can 
write this last expression as 

max{ [logp aj, [log p fej} 

a + b= 22 P k d p 2 ( 0 , d p (k, a) + d p (k, b)) ( 12 ) 

k=— oo 

Therefore, by using Eq. ([8]) 

max{ Llog p aj, [logp 6J} max{ [log p aj, [log p b\ } 

a+b = 22 P k d p (0,d p (k,a) + d p (k,b))+ 22 P k+l d p (1, d p (k, a) + d p (k, b)) 

k=—oo k=— oo 

(13) 

Thus, we see that the ordinary sum of two real numbers can be split into the bitwise addition modulo 
p of a and b 

max{[log p aJ,Llog p b|} 

a+ p b= 22 p k d p (0,dp(k,a) + d p (k,b)) (14) 

k=— oo 

plus a term which provides the total contribution of the carries 

max{ [logp aj, [logp 6J} 

a+pb= 22 P k+1 d p (l,d p (k,a) + d p (k, b)) (15) 

k=— oo 

Note that if the kth-digits of a and b are such that d p (k, a) + d p (k, b) < p there is no contribution to 
the carries coming from the fc-th digit since in that case d p (1, d p (k, a) + d p (k, b)) = 0. We thus have 
proved 

a + b = (a +p b ) + (a+pfr) (16) 

Our goal in this paper is to find all possible bitwise operators which, as is the case of + p , act on 
two arbitrary real numbers. This goal is attained by first constructing all discrete operators that are 
laws of composition (magmas) of S x S -» S. We denote such operators by 2 -R p (xo> x\). They act on 
integers xq and x\. both in S. The operator can be defined through a table as 
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2-Rp 

0 

1 

...k... 

p — 2 

P~ 1 

0 

ao 

ai 

... a,k • ■ ■ 

dp— 2 

dp—i 

1 

dp 

&P+1 

• • • ttp+k • • • 

a 2p ~ 2 

a 2p ~ 1 

2 

a 2p 

02p+l 

. . . d2p-\-k • • • 

d3p- 2 

«3p-l 

...k... 

^kp 

Q-kp+l 

• • • &kp+k • • • 

^(k-\-l)p—2 

1? 

+ 

•S' 

1 

p — 2 

a p(p- 2) 

a p(p— 2)+l • 

. . 2)-\-k • 

• ^p(p— 1)—2 

^ J p(p—1)—1 

p — 1 

a p(p- 1) 

a p{p~ 1)+1 ’ 

. . dp^p_ • 

a p 2_2 

tip 2 — 1 


where xo is listed in the columns and x± in the rows, and the result a xo+pxi 6 S, is given by the table. 
Clearly, we have for any a ne pj that n 6 [0 ,p 2 — 1]. Since all a n ’s are integers £ [0 ,p — 1] this table 
provides the most general description of a law of composition (also called a magma) [571 [38] . The 
integer number R entering in the label of the operator 2 R p is defined as 


p 2 - 1 

R=^2 a nP n ■ 

n =0 

2 

Thus, R £ [0,pP — 1]. Therefore, it is clear that, from the above table 

/ p 2 -i \ 

2 R p (x 0 ,x 1 ) = a Xo+ p Xl = dp I xo +pxi, 22 p n a n = d p (x 0 + px±,R) 

\ n =0 / 


(17) 


(18) 


since the (xq +pxi)-th digit of the number Yln=o P na n a x 0 +px 1 - This last expression in terms of 
the digit function most concisely expresses any possible magma. Let now u and v be arbitrary real 
numbers. The generalized bitwise operator b g ( 2 R p ',u,v), with q > 2 acting on these real numbers is 
defined as 


b q ( 2 R p ,u,v) = 22 <l k 2-ftp (dp(k,u),d p (k,v)) 

k =—00 
kmax 

= 22 qkd P ( d p(k,u) + pd p (k,v),R) 


(19) 


k =—00 


where k ma x = max{ [log p rij, [log p uj}. This is the expression that we shall analyze in the rest of the 
article. 

Example: Let us find what is the numerical value of b 2 ( 2142 ; 5,11). We have that 5 and 11 are given 
respectively in base p = 2 as ’0101’ and T011’. Now the operator 2142 has vector (ao, ai, a 2 , < 23 ) = 
( 0 , 1 , 1 , 1 ) which means that it returns one when the digit k of u or v or both is one and zero otherwise. 
Then, the resulting string is (1 1 1 1) and, therefore, b g ( 2 M 2 ; 5,11) = 1 • q° + 1 • q 1 + 1 • q 2 + 1 • q 3 . 
Thus b 2 ( 2142 ; 5,11) = 15. 

Let us now also see how the operator + p is a particular instance of Eq. (19). We note that if we 
take a n = d p ( 0 , d p ( 0 , n) + d p ( 1 , n)) so that 


p 2 -1 


r = 22 pHd p (°> d p (°> n ) + d p ( l > n )) 


( 20 ) 


n =0 
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we then have, from Eq. (18) 


p 2 -1 


d p (d p (k,u) + pd p (k,v),R) = d p d p (k,u) + pd p (k,v), ^ p n d p (0,d p (0,n) + d p (1 , n)) 


n =0 


= dp (0, dp (0, d p (k, u) + pd p (k, v)) + d p (1, d p (k, u ) + pd p (k, v))) = d p (0, d p {k, u) + pd p {k, v)) 


Therefore, by taking q = p in Eq. (19) and replacing this expression we obtain 

kmax kmax 


b p ( 2 R p ,u,v) = P k d p (d p (k,u) + pd p (k,v),R) = ^ p k d p (0, d p (k, u) + pd p (k, v)) 

k =—oo 

= u+ p v 


k =—oo 


( 21 ) 


as can be seen from Eq. (14). Thus, the addition modulo p of two real numbers u, v is a particular 
instance of a generalized bitwise arithmetic operator h q (2 R p ',u,v) with q = p and with the magma 
2 R p having code R given by Eq. (20). 

We note that our expression Eq. (18) for any magma in terms of the digit function can be gener¬ 
alized to an arbitrary Wary operator jv-Rp : S N —> S as 


r iV-l 


P N - 1 


fN-l 


nR P (xq , • ■ -,XN-l) = d P E pk ' Xk - E p71an = d P E pkxkl R 


( 22 ) 


k =0 


n =0 
n N 


,k =0 


N _ 1 Jy 

where now R = Yln=o a n'P n is an integer E [0,^" — 1]. This operator can be used to define a 
generalized bitwise operator acting on N real numbers uq, u±, ..., u^-i as 

max{ Llog p u 0 J, Llog p Ujv- 1J } 

h q ( N R p ,u 0 ,...,u N -i) = ^2 N Rp(d p (k,u 0 ),...,d p (k,u N -i))q k 

k =—00 

max{ [log p u 0 J, Llog p Mjv- lJ } 


m -1 


E 

k =—00 


''/'dp p n d p (k,u n ),R 


(23) 


v n =0 


This function acts in the following way: it takes the digit k of each number uq,. . rtjv -1 in radix p 
and introduces the N-ary operator nR p given by Eq. (22) on the N digits. Then, the resulting string 
of digits is interpreted as a number in radix q (if > p) or simply, in general, as a polynomial in powers 
of q. Thus, b g is a non-negative real number. 

We now prove one of the main results of this article. 

Theorem 2.1 (Self-affinity of the generalized bitwise function.) The function b q defined 
by Eq. (23) has a self-affine property, so that 

bq ( nR p ]u 0 , • • .,u N - 1 ) = p~ hlq b q ( N R p ;pu 0 ,... ,pu N - 1 ) (24) 

so that H = log p q defines a ’roughness exponent’ of a hypersurface in N dimensions, cf. Eq. |7p. 









Proof: We have 


max{ [log p pu 0 J,., Llog p pu N _ 1 J } 

bg ( N Rp-,pu 0 ,... ,pu N -i) = ^2 NRp(d p (k,pu 0 ),...,dp(k,pu N -i))q k 

k=— oo 

1 +max{ L lo §p “0 J, • • ■ > U°gp “ N- 1J } 

= ^2 N Rp(d p (k-l,u 0 ),...,d p (k-l,u N ^ 1 ))q k 

k=— oo 

max{ [log p -uojL lo gp «JV-1J} 

= ^2 NRp{d p {k',u 0 ),...,d p (k',u N -i))q k+1 

k'=— oo 

max{ [log p ttojL lo gp “JV -1J} 

= p lo s P i ^2 N R p (d p {k',u 0 ),...,d p (k',u N -i))q k ' 

k'=— oo 

= p log p q h q ( n R Pi uq,...,u N -i) ( 25 ) 

where we have used the scaling property, Eq. 0 in noting that 


d p (k,px) = d p (k — l,x) 


(26) 


This proves the result and that log p q = H is the roughness exponent of the self-affine surface by 
comparing this latter expression with Eq. 0 . □ 

The crucial property established in general by the above theorem proves that the operator x + p y 
in Fig. [I] has indeed a self-affine property and deserves to be called a fractal in its own right. 

We also have the following result. 

Theorem 2.2 (Coarse-graining theorem). Let nRp be a N-variable discrete operator 
in radix p acting bitwise on the digits of the reals uq, u\, ..., un-i- Then, the function 
b q (R p ;uo, u\, ..., un-i) defined in Eq. (23) is smoothened as q is increased. Furthermore, if 


q » p the following expression is asymptotically satisfied 


bg ~ 


b T 


ph 


where k max = max{ [log p u 0 J, L lo gp «ij, • • ■, |>g P “JV-iJ }• 

Proof: Since h q is non-negative by definition, from Eq. (Jd]) we have h q = J2k=-oo ( 1 
coefficients accompanying powers of q in Eq. 


k d q (k,b q 


(27) 


All 


(19) are non-negative integers. 


Thus by increases 

monotonically with q. Since all coefficients are bounded (and all are independent of q) and their 
dependence on the coordinates u n , n = 0,..., N — 1 is governed by k it is clear then that, as q increases, 
terms with k larger become more significant and the sum displays a longer spatial variability. As a 
consequence of this, details are gradually lost and coarser surfaces are obtained. Now, note also that 
b p = Ylk'= a -oo P k d p (k, bp). When q » p is large only the most significant digit of b g becomes relevant 
and we have 


by ~ 


c H ( k 


max^ bg) 


Now, since from Eq. (23), for q > p we have 

d q (k,b q ) = d p (k, bp) = d p ^ p n d p (k,u n ),R 


/N -1 


(28) 


(29) 


vn=0 
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FIG. 3 1 Two views of the function b 3 (2 139033; u, v ) in a portion of the plane with u £ [0, 200] and v £ [0, 200] 


and the digits of b g are the same as those of b p and, hence, belong to the set S. Then, from Eq. 
(28) and, by using the definition of the digit function, Eq. ^ and noting that \h p /p kmax+l \ = 0, we 
finally obtain 


b q ~ 


q kmax d q (k max ,b q ) = q kmax d p {k max ,b p ) = 


bp 

P^max 


(30) 


which proves the result. □ 

The coarse graining theorem implies that b g asymptotically behaves like the most significant digit 
of bp for q » p, after a trivial rescaling (i.e. dividing by q kmax y Since all information contained 
in the less significant digits of b p is lost in b q , the latter exhibits no fine details compared to the 
former and if one takes the numbers uq, u\, ..., u^-i as ’coordinates’ only the broadest patches in 
the space spanned along these coordinates survive: Only the operations on the most significant digit 
of the numbers uq, ui, ..., un~i become relevant. 


III. PATCHWORK QUILTS: SELF-AFFINE SURFACES AND THE ROUGHNESS 
EXPONENT 


To better substantiate these results and study particular examples, it is better to fix N = 2 in Eq. 


(23) i.e., to work with Eq. (19), so that we can represent such operators as functions taking values on 
the plane given by coordinates u and v. 
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FIG. 4: From top to bottom and left to right, the functions b 2 ( 262 ; u, v), b 3 ( 274173 ; u, v), b 3 ( 294073 ; u, u) and 
b 5 ( 2 13427417 5 ; u, v), with u S [0,100] and v S [0,100] 


In Figs. [ 3 ] and is plotted as a function of (u, v ) and the operator indicated in each case. 
Patchwork quilts are obtained in every case when projecting the 3D object on the plane (the figures 
just only show a finite portion of the plane). Such representations on the plane are useful to get a 
glimpse of the behavior of the operator 2 R P within b g . For example that 262 is commutative can be 
clearly seen from the reflection axis on the southwest-northeast diagonal exhibited by b 2 ( 262511 , v) 
(see Fig. [4j top left). The different patches give information on how h q organizes in the plane. The 
tridimensional aspect of these surfaces, an example is shown in Fig. [4j reminds the towers of the 
’exotic’ fractal surfaces reported by Indekeu et al. obtained through random deposition [20] . 

The roughness exponent obtained above has some of the nice properties reported in experimental 
studies on fractal self-affine surfaces [29j [33] [35]: A larger value of the exponent corresponds to a 
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FIG. 5: Gradual spatial coarsening of the bitwise operation function h q as the parameter q is increased. Shown is the function 
b q (2139033; u, v) for all integer values of q E [3,11] as indicated besides each panel (u G [0,100] and v E [0,100]). 


coarser (i.e. ’smoother’) surface. This is a result of the coarse graining theorem in the previous section, 
which we illustrate in Fig. [5] for the case q > p (the case q < p follows a similar trend although not 
that markedly). In Fig. [5j the function h q (13903s; u, v ) is plotted for the values of q e [3,11] (p = 3) 
shown in the Figure. We observe that for q increasingly large, details are gradually lost and only the 
coarser structures survive. In the Figure it is observed that for q = 10 the asymptotic limit is already 
being approached. Remarkably, for q prime details sometimes reemerge (but in a coarser scale). This 
rather subtle phenomenon poses an interesting number theoretical problem for which we have not yet 
a solution, and which shall be addressed elsewhere. 


Increasing q leads to coarser structures. If q is not asymptotically large the structures at different 
scales overlap creating complex patterns that are coarser on the average. We can then ask the question 
whether is it possible to sharply discriminate the coarser structures without significantly changing the 
value of the function (beyond, of course, the change induced by neglecting the fine details). This is 


12 




























4 

■ 1 

II 

. B=B: 

• B 

: 


■LLP 1 ■ ■ ■ 

m m ■ i ■ 


wist;::*; 

i I ? » B:B:. E 
ru i 

I ■ • • B:B:1 ii 

k:b:sk:b:s 
l B1B:« B!B:S b 




■ 


I I 

■ i 

■ I 

:: 

■ 


■ . ■ 
I ■ I 
III 
I . I 
I - I 
III 
I . I 
I 


FIG. 6 ; Action of the coarse graining operator on a self-affine surface, as described by Eq. (|3lj). Shown is the quantity 
3~ D [3 D h3 (298153; u, u)J for the values D = 0 (left), D = —1 (center) and D = —2 (right) with u E [0,100] and v E [0,100]. 


accomplished by the coarse graining operator acting on b g , i.e. q u Yq u b q { 2 Rp]u,v)\. By applying 
this operator to Eq. (19) we obtain 


q D [q D h q ( 2 R P ;u,v)\ = ^ q k d p (d p (k,u) + pd p (k,v), R) 


(31) 


k=—D 


The error committed in approximating b 9 by such truncation is thus strictly smaller than q~ D . Since 
each term accompanying a power k of q contributes a spatial structure (that is coarser for k larger), 
making D more negative selects the coarser structures while making the value of q~ D [q D h q (2 R p ; u, u)J 
to slightly decrease. 


The action of the coarse graining operator, as described by Eq. (31), is shown in Fig. [h] for the 
function 3~ D [3 D h3 ( 298 I 53 ; u, u)J for the values D = 0 (left), D = — 1 (center) and D = — 2 (right). 
We see that the coarse graining operator merely selects broader patches, eliminating the fine details. 

The structure of the coarse graining operator is also interesting in the sense that it reveals that the 
operation of multipliying by an scalar, i.e. q~ D , and taking the floor operator [• ■ -J do not generally 
commute: q~ D • -J 7^ \_q~ D • • -J • This lack of commutativity is behind the loss of information implied 
by the coarse graining. In [32] we have discussed the implications that, in our view, this fact has for 
the foundations of quantum mechanics. 


IV. CONCLUSIONS 

In this article a new general method to mathematically design fractal surfaces with self-affine 
properties has been presented. We have discovered these structures hidden within operations as simple 
as the ordinary sum of real numbers. Thus, we speculate that such structures may be of general interest 
in the mathematical modeling of complex physical systems. Our approach is based on digital calculus, 
which we have introduced to describe rule-based dynamical systems, as cellular automata m and 
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substitution systems |13j . This approach is also the cornerstone of a recent formulation of quantum 
mechanics [32] and of our method to derive fractal decompositions of conserved physical quantities 
[33] . We have rigorously proved that bitwise operators acting on real scalar fields on each point in 
the plane M 2 lead to surfaces with self-affine properties. We have obtained the roughness exponent for 
these surfaces and shown that it has the characteristic nice properties expected for such an exponent, 
as widely reported in experimental systems: The larger this exponent, the coarser the resulting surface 


The fractal surfaces here obtained come under the general name of ’patchwork quilts’ because we 
are following the acknowledgment made by Donald Knuth in [39] to one pattern designed by D. Sleator 
in 1976 (unpublished result, displayed on p. 136 of [39]) resembling the ones obtained here. However, 
no relationship of such pattern to fractals and self-affinity seems to have ever been attempted and no 
systematic construction of such patterns is found in the literature. In this article, Eq. (23) constitutes 
a universal form embodying all conceivable bitwise operators in any radices p and q of standard 
numeral systems and any number of real scalars N, the different digits representing independent bits 
of information at different scales. 

We speculate that our surfaces can be related to free energy landscapes of clusters and bionrolecules 
EQ! obtained under the basin hopping scheme m- An analysis of this connection shall be investigated 
elsewhere. 
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